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ABSTRACT Robertson's Mutator {Mu) systenn has been used in large scale mutagenesis in maize, exploiting KEYWORDS 
its high mutation frequency, controllability, preferential insertion in genes, and independence of donor Mutator 
location. Eight Mutator elements have been fully characterized (Mu7, Mu2/Mu1.7, Mu3, Mu4, Mu5, Mu6/7, transposon 
Mu8, MuDR), and three are defined by TIR {Mu10, Mu11 and Mu12). The genome sequencing revealed Mutator13 
a complex family of Mu-like-elements (MULEs) in the B73 genome. In this article, we report the identification mutagenesis 
of a new Mu element, named Mu13. Mu13 showed typical Mu characteristics by having a ~220 bp TIR, maize 
creating a 9 bp target site duplication upon insertion, yet the internal sequence is completely different from transposon 
previously identified Mu elements. Mu13 is not present in the B73 genome or a Zea mays subsp. parviglu- tagging 
mis accession, but in W22 and several inbreds that found the Robertson's Mutator line. Analysis of mutants cloning 
isolated from the UniformMu mutagenic population indicated that the Mu73 element is active in trans- 
position. Two novel insertions were found in expressed genes. To test other unknown Mu elements, we 
selected six new Mu elements from the B73 genome. Southern analysis indicated that most of these 
elements were present in the UniformMu lines. From these results, we conclude that Mu13 is a new and 
active Mu element that significantly contributed to the mutagenesis in the UniformMu population. The 
Robertson's Mutator line may harbor other unknown active Mu elements. 



Mutator (Mu) transposable elements are a major class of class II trans- 
posons identified in maize by Donald Robertson (1978, 1981). The 
two-component system, one autonomous MuDR and many nonauton- 
omous Mu elements, was exploited for efficient mutagenesis in maize. 
High copies of the elements offer a high forward mutation frequency, 
whereas limited copies of MuDR allowed turning off the transposition 
by removing the element through segregation (McCarty et al. 2005). 
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Preferential transposition into gene rich regions by Mu elements en- 
hances mutagenesis frequency. And transposition not limited to linked 
loci facilitates genome wide mutagenesis. For these reasons, several 
mutant populations in maize were created by using the Mu system 
(Bensen et al 1995; May et al 2003; Raizada 2003; McCarty et al 2005). 

The well-characterized Mu elements {Mul to Mu9/MuDR) were 
discovered exclusively in maize. Subsequent molecular analyses and 
genome sequencing revealed that Mu elements are present in plants 
(Lisch 2002), fiingi (Chalvet et al 2003), bacteria (Eisen et al 1994), 
protozoans (Pritham et al 2005), and metazoans (Hua-Van and Capy 
2008). Based on sequence similarity, these elements are classified as 
Mu-like elements (MULEs). MULEs belong to a superfamily of trans- 
posons with complex members and diverse sequences. Typical char- 
acteristics of this family include a conserved 50-200bp terminal 
inverted repeat (TIR), unrelated internal sequences between the TIRs, 
and creating a 9bp target site duplication (TSD). In contrast, all the 
previously identified Mu elements from maize (Mul-Mu9/MuDR) 
carry a ~220bp TIR that is highly conserved. Transposition activity 
of the elements is thought to be associated with the TIR sequences. 
Inactive elements carry mutated TIRs. 

Different from Ac/Ds and Spm/dSpm transposable elements 
where the non- autonomous elements are deletion derivatives of 
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the autonomous elements, the internal sequences between TIRs 
among Mu elements are often unrelated. Some Mu internal sequen- 
ces showed high similarity to host genome, suggesting a possible 
gene capture in the formation of these elements. This class of Mu 
elements was classified as Pack-MULEs (Jiang et al 2004). About 
262 Pack-MULEs were identified in the B73 genome (Schnable et al. 
2009). Because promoters are found in the TIRs, Mu internal 
sequences may be transcribed in convergent orientations (Hershberger 
et al. 1995; Lisch 2002). Hence, it was suspected that some of the 
Pack-MULEs may have regulatory fianction, as antisense transcripts 
may interfere with expression of the endogenous genes (Lisch 2005; 
Juretic et al. 2005). 

Transposition of all Mu elements required the presence of an 
active MuDR element. The MuDR element contains two genes, mudrA 
encoding a transposase (MURA) and mudrB whose product (MURE) 
is of unknown function. MURA showed high similarity to bacterial 
transposase and the virus integrase (Walbot and Rudenko 2002); 
hence, it is essential for transposition. Transposable elements con- 
taining only mudrA-like genes were found in both monocots and 
eudicots (Saccaro et al. 2007). The mudrB gene is only present in the 
genus Zea (Lisch et al. 2001). Jittery, an autonomous transposon 
identified in maize, contains a mudrA-like gene, but with TIR 
sequences distinct from Mu elements (Xu et al. 2004). Jittery 
exhibited high frequency of excision, causing somatic and germinal 
reversion, but apparently lost its activity for new insertions. Trans- 
position of Mu elements employs two distinct mechanisms. In so- 
matic cells, transposition mostly uses a "cut-and-paste" mechanism. 
The element cuts itself and reinserts it in a new locus elsewhere in 
the genome. High-frequency excision of Mu elements is restricted to 
the late stage of cells in development during organogenesis. In ger- 
minal cells, Mu transposition uses a "replicate-and-insert" mecha- 
nism where the element replicates just before meiosis or in the 
gametophytes and inserts in a new locus in the genome. Conse- 
quently, "cut and paste" transposition does not increase the copy 
number, whereas "duplicate-and-insert" transposition does. Excision 
of a Mu element left a footprint of the 9bp TSD, which sometimes 
restored the function of the donor gene such as in hzl-mum9 
(McCarty et al. 2005). 

Prior to the sequencing of the maize genome, eleven Mu elements 
were reported in maize, of which eight were characterized by fiall sequen- 
ces, (/.e., Muh Mu2 IMul.7, Mu3, Mu4, Mu5, Mu6/7, Mu8 and MuDR) 
(Bennetzen et al. 1984; Taylor and Walbot 1987; Talbert et al. 1989; 
Fleenor et al. 1990), and three were indicated by TIRs (Dietrich et al. 
2002). The sequencing of the B73 genome revealed a surprisingly com- 
plex view of the Mutator family, which accounts for approximately 1% of 
the 2.3 gbp genome (Schnable et al. 2009). These include MULEs, Pack- 
MULEs, and SOLOs that contain only one TIR. Many of these elements 
contain a shorter TIR, suggesting that these elements may have lost the 
capacity for transposition. In this study, we report a new Mu element, 
Mul3y which was identified from the UniformMu population, a deriva- 
tive Mu active line from the Robertson's Mu line. Mul3 exhibits typical 
Mu characteristics and is active in transposition. It contributes signifi- 
cantly to mutagenesis. The finding of Mul3 adds to the active Mu 
reservoir and facilitates cloning of causative insertions in the Mu tagged 
mutants in phenotype-driven forward genetics in maize. 

MATERIALS AND METHODS 
Genetic stocks 

The maize lines used in this study were derived from the UniformMu 
population, a Mutator line with the mutable hzl-mum9 anthocyanin 



biosynthetic gene introgressed into the W22 genetic background 
(McCarty et al. 2005). The teosinte lines Zea mays subsp. parviglumis 
(Accession: PI 384061) and Zea mays subsp. mexicana (Accession: PI 
566684) were provided by the Maize Genetic Stock Center. Other 
inbred lines (W22, B73, Mol7, M14, Q66, Q67, B77, and B79) were 
generously provided by Donald R. McCarty (University of Florida). 

Cloning of Mu13 from UniformMu population 

The Mul3 transposable element was amplified by a pair of primers 
(5'-CTGCTCCTGTGCTATCCTCC-3' and 5'-ACCAAACCAACA 
AGAGCCTG-3') flanking a Mul3 insertion in a gene coding for 
a putative plastid Sigma factor 3 (ZmSig3). Template DNA was isolated 
from line 03S-4081-01, homozygous for the insertion. As Mu elements 
carry a long terminal inverted repeat (~220bp), it interferes with PGR 
amplification. We tested different conditions with DNA polymerases 
of various sources. ExTaq (TaKaRa, Japan) and ThermalAce DNA 
polymerases (Invitrogen, USA) yielded successfiil amplification. The 
PGR reaction was composed of 20mM Tris-HCl pH 8.4, 50 mM KCl, 
2mM MgCl2, 200 |jlm of each dNTP, 100 nM each primer, 5% DMSO, 
and 1 U of DNA polymerase. PGR conditions were 96°G/3min for 
initial denaturation, 8 cycles (95°G/30 sec, 62°G/30 sec, 72°G/2min) 
followed by 30 cycles (95°G/30 sec, 58°G/30sec, 72°G/2min), with final 
extension at 72°/10min. The PGR fragment was purified from gel by 
gel extraction kit (Zymo Research, USA), ligated into pGR4-TOPO 
(Invitrogen, USA), and sequenced. 

Selection and cloning of new Mu elements in B73: A conserved 
200bp Mu TIR sequence based on known Mu elements (Mul to Mu9/ 
MuDR) was used in a BLAST search of the GenBank maize sequences, 
with a cut-off E value of < e~^^. Within this collection, the known Mu 
elements were identified by a BLAST search with the internal se- 
quence of each Mu element. Identical sequences were clustered using 
BLASTGLUST (ftp://ftp.ncbi.nih.gov/blast). The resulting collection 
was analyzed for lefl;- and right-TIR in terms of orientation and ho- 
mology, as well as the presence/absence of a 9bp host sequence direct 
duplication. 

We amplified the internal sequences of six new Mu elements that 
carry highly conserved TIR at both ends. The primers were listed in 
Table I, and the PGR conditions were similar to those present in the 
amplification of Mul3. The internal sequences were cloned in pGR4- 
TOPO and sequenced. 

Selection of UniformMu mutants for Southern blot analyses: 

UniformMu mutants segregating for visible mutant phenotype of 
embryo defective (emb), small kernel (smk), empty pericarp (emp), 
shrunken (sh), and defective kernel (dek) were randomly chosen. The 
18 mutants were 06S-6001 (smk), 06S-6002 (emp); 06S-6004 (de- 
fective kernel, dek); 06S-6005 (emp); 06S-6016 (smk); 06S-6018 (dek); 
06S-6019 (smk); 06S-6020 (emb); 06S-6023 (emp); 06S-6026 (smk/ 
dek); 06S-6029 (smk); 06S-6032 (smk); 06S-6033 (emb); 06S-6034 
(dek); 06S-6044 (dek); 06S-6045 (emp); and 06S-6055 (sh/smk). Each 
DNA was extracted from seedlings of three individual ears that were 
genotyped based on the seed phenotype. All these ears did not exhibit 
active MuDR activity, as indicated by the mutable bzl-mum9 antho- 
cyanin biosynthetic marker. AU these lines were back- crossed with 
W22 twice. 

DNA extraction and Southern analysis: Genomic DNA was isolated 
by a urea-phenol-chloroform-based method. Ig fresh weight of leaf 
tissues was ground in liquid nitrogen and extracted with 5 ml of DNA 
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extraction buffer (7 M urea, 0.3 M NaCl, 50 mM Tris-HCl, 24 mM 
EDTA, and 1% sarkosine, pH 8.0). After mixing with 4 ml phenol- 
chloroform-isoamyl alcohol (25:24:1), the extraction was carried out 
with gentle shaking for 30 min at room temperature. The mixture was 
separated by centriftigation at 4800 x g for 15 min. The aqueous phase 
was transferred to a new tube and mixed with 0.1 volume of 3M 
sodium acetate (pH 5.2) and 3.8 ml isopropanol. DNA was pelleted 
at 4800 X g for 5 min, washed with 70% ethanol, and dissolved in TE 
buffer (10 mM Tris-HCl, 1 mM EDTA, pH 8). Approximately 10 |jLg 
genomic DNA was digested with appropriate restriction enzymes at 
37°C for 6 hr. The DNA was resolved on a 0.7% agarose gel, dena- 
tured, and blotted onto a Hybond-N membrane (GE Healthcare). The 
membrane was cross-linked and hybridized. The probe was labeled 
with Ready-To-Go DNA labeling beads and purified with ProbeQuant 
G-50 micro column (GE Healthcare). 

The probes used for Southern analyses for Mul4-Mul9 were am- 
plified from the B73 genome by PGR with primers listed in Table 1. 
The primer anchor positions with respect to TIR and probe sequences 
are listed in supporting information. File SI. For MulS, it was derived 
from UniformMu by PGR with a single primer (5'-ATGAATGTGGT 
GTGAGGGTTTAGGGT-3') that was anchored in the TIR region. 

Bioinformatics analysis: Sequence alignments were carried out using 
the GLUSTALW algorithm available online (http://workbench.sdsc. 
edu/). For phylogenetic tree construction, the phylogenetic tree files 
from GLUSTALW analysis were imported into a Tree View program 
(http :/ / taxonomy.zoology.gla.ac.uk/rod/ treeview.html) . 

RESULTS 

Identification of Mu13 element 

In a large scale extraction of Mu flanking sequences from mutants 
isolated from the UniformMu population (McGarty et al. 2005), a Mu 
element was found inserted in a gene coding for a putative plastid 
Sigma factor 3 (ZmSigS). The element was inserted in the third exon 
of the gene (refer to Figure 4, Accession no. GG893004). We cloned 
the Mu element and found that the element is 1494bp long, contain- 
ing a 223bp TIR with an 88% identity to the consensus of previously 
known Mu TIRs (Figure 1). The left and right TIR showed a higher 
identity (92%), which is a general feature of the Mu elements. A 
9bp direct target site duplication was found at the insertion site. 
The internal sequence of this element is completely different from 
any previously identified Mu elements (Bennetzen 1984; Taylor & 
Walbot 1987; Ghen et al 1987; Talbert et al 1989; Fleenor et al 



^ Table 1 Primers used in amplification of selected Mu elements 
in B73 



Primer 




Sequence 


Mu14-iF1 


5' 


-CTCTTCCCCACACCTATTGC-3' 


Mu14-iR1 


5' 


-GAGATGCTCCGCGATTACAT-3' 


Mu15-iF1 


5' 


-TAAGGTGATTTGCTCGGGTC-3' 


Mu15-iR1 


5' 


-TCTCTTGCTTCTCCGTCTCC-3 ' 


Mu16-iF1 


5' 


-CACCGTCAGGCTTAACAACA-3' 


Mu16-iR1 


5' 


-CGGTGAGTTCTCCTCCTCTG-3' 


Mu17-iF1 


5' 


-CTCAGCGAACTCTGGCACAC-3' 


Mu17-iR1 


5' 


-CACTCCTCTCCGTCTCCG AT-3 ' 


Mu18-iF1 


5' 


-TTGGAGGTGTCGGTAGTGAGC-3' 


Mu18-iR1 


5' 


-ACAGCTCTTGCGTCTCCTCTG-3 ' 


Mu19-iF1 


5' 


-ATTGGAGTGCTCTCGGGGT-3' 


Mu19-iR1 


5'- 


-AG AGCTCGGTCTCAGGCATTA-3 ' 



1990; Hershberger et al 1991). Searching the GenBank and the nearly 
completed B73 genome did not find the presence of this element. In 
light of the partially characterized Mu elements, MulO, Mull, and 
Mul2 (Dietrich et al 2002), we designated this element as Mul3 
(Accession: HQ698272). 

Bioinformatic analysis revealed that the Mul3 internal sequences 
contain two open reading frames (ORF). The conceptually translated 
protein sequences of these two ORFs showed high similarity to a maize 
protein that was annotated as nucleotide binding protein (accession: 
AGG25371, GRMZM2G3 17614). Further analysis revealed that it enc- 
odes a WD40 protein, containing seven WD repeats. As indicated in 
Figure ID, the first highly similar ORF started from the first methi- 
onine and covered 73 amino acid (aa) residues in length. This region 
shared an 88% identity with the maize WD40 protein, and a similar 
identity with apparent orthologs in sorghum (Sb01g008680) and rice 
(Os03g0738700, also identified as Os03g52870, annotated as trans- 
ducin family protein). The second highly similar region (95% identity) 
was about 42 aa long and coincided with the first repeat of the WD40 
protein. In the maize WD40 protein, these two regions were separated 
by 26 amino acid residues, which were not found in the Mul3. This 
maize WD40 gene was expressed as indicated by ESTs, suggesting that 
it may be a functional gene. Another WD40 gene on maize chromo- 
some 5 (GRMZM5G852097) is apparently a syntentic paralogous du- 
plicate of GRMZM2G3 17614, which is also probably functional. 

Mu13 is active in the UniformMu population 

The insertion of MulS in a fianctional gene in the UniformMu pop- 
ulation suggested that it was active in transposition. This insertion was 
not present in the parental lines that gave rise to the mutant. It is 
known that Mu elements are not equally active. Mu4, Mu5, and Mu7 
were less active than the other known ones (Talbert et al 1989), and 
so far most genes cloned by transposon tagging were inserted by Mull 
2, Mu3, Mu8, and MuDR. To assess the MulS transposition activity, 
we analyzed 18 UniformMu seed phenotype mutants randomly se- 
lected from a large set of available seed phenotype mutants. For each 
mutant, seeds showing no MuDr activity (lack of somatic transposi- 
tion indicated by the hzl-mum9 marker gene) were selfed to produce 
an F2 mutant segregating family. The genotype of each F2 individual 
was scored by examining the ear. DNAs from three F2 individuals of 
either wild type (not segregating mutant phenotypes, N) or segregat- 
ing seed mutant phenotype (S) were pooled separately and analyzed 
by Southern hybridization. As shown in Figure 2, hybridization with 
a Mul3 probe detected seven new MulS insertions, as indicated by the 
appearance of new MulS containing fragments. Because the MulS 
probe used in this analysis contained 80bp TIR sequences, it cross- 
hybridized with related Mu elements and produced weak signals. The 
MulS signals were strong. Three MulS containing fragments (4.5kb, 
5.8kb, >12kb) showed uniform presence in all the members, suggest- 
ing that they are apparently parental. When the same blot was hy- 
bridized with a MullMu2 specific probe, comparable numbers of 
Mul/Mu2 insertions were detected (Figure 2B). Some of these inser- 
tions were unique to individual lines, suggesting new transposition by 
Mul/Mu2 as well. This result indicated that MulS is active in the 
UniformMu population. 

Mu13 presence in maize inbred lines and teosinte 

The UniformMu mutagenic population was derived from introgress- 
ing Robertson's Mw-active line into W22 genetic background 
(McGarty et al 2005). Hence, the MulS element can be derived from 
either W22 or Robertson's Mw-active line. To determine the presence 
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GAGATAATTGTCATTATGGAAGAAGAGAGGAGGGGATTCGACGAAATGGAGGCGTTGGCGTTGGCTTCTCTGTTTTAGAGACGCACGCGACA 
GACAAACTCTAAAATGGATACGAGACAACACTTAGGGCTACGTAAACCGGTATCAATGTCCTGTCACCGTTTACCGTTCACGTGCACAGACG 
CCGTCTGCGTACTCTTCTCGTCTCCGTCGCTTCTCCTTC ATCTCGGAATTGGAGTGCCCTCGCCCTCGCCGTCGCTGTCCTGTCATCTGTTG 
TTCCCCATGTCTAAGTCCGAGATGAACTGTACCGGAGCTCCACCACGTCCTCCGGCAAGAAACATATGAAAGACAGCGACGGCGAGGGCGCC 

MKDSDGEGA 
GGCGGTGGCGTCAACCGCTCGCACCCGTCCAACAACCTTTCACTGCCGGCTCCATACTCCGACCCGAACCTCCAGTTCTCGGGTGGGACGGA 
GGGVNRSHPSNNLSLPAPYSDPNLQFSGGTD 
CGATGAGTCCTCGAACCGGGACAACAGCTCCTCTGCCACCGGCGGCGCGAGCCCGAGGTACTACTGCGACTACCCGTCCAGCTTCAGCGGCG 
DESSNRDNSSSATGGASPRYYCDYPSSFSGE 
AGTGCTCCTCTGCCACCGGCGGCGCGCAGGAAGGAAAAATAAAAAAAAAATCTACTGCAGCCGAGCCGCCGCCCCTTGGTGAGATTTTTTCC 
C S 

TGCGCAACGTGTCCCTGCTCGCTTCTTCTAGAGCTCGAGGTCATCTTCTCCAGCCTCCATCTCCCTCACAGCCGCGCGCTGCTCCAGCTTCT 
TCCCCTGCGCGTCGCGCACTTCTTCCCTGCTGGCGCCGCGCCCTGAGCTTCTTCGCCACCGCGCCCCTGGTTCCGGCCCCGGCTCGACGGCC 
TTCTTCCCAGGCGAGCTCCCGGCGTTGGATCCGCTGCGAGCTCCCCTGCAGCTCGTCCCCATCTTTCTGGCCGCGCTGAAGTCCCGCTTCGT 
GGGAGTCTGGAGCAGCCGTGGACGCGCCCGCGCTGTCCCTGCTCGTGTCTTCTCCCATGGCGCCTGGCTAGATCTCCCTGCTGCGCTCGAGA 
TAATTGCCATTATGTAAGAAGAGAGAAGGGGATTTGACGGCCCGCCCACGATGGCGCCGGAGACGAGCCTCATCAGCTCGCTGGTGATGGAG 

GPPTMAPETSLISSLVME 
GAAGGCCACATCTACTCGCTCGCCGCCAAGGGCGACACGCTGTACACCGGCTCGGACAGCAAGAACATCCGAGCCGGGTATCTCGCCGGCGT 
EGHIYSLAAKGDTLYTGSDSKNIR 

GCGCTCCTGCGCCGTGGGGTCGCAAGACGGAGGCGTCTCTCATGCGCTCCGCACGGGCCTCCTCGTCGCCGGCGAAGGAGAAGCGACGGAGA 
CGAGAAGAGTATGCCAGACGGCGTCTGCCCACGGGAACGGTAAACAGTGACAGGACATTGATACCGGTTTATGTAGCCCTAAGGGCTGTCTT 
GTAT CCATTTTGGAGTTTTGTCT GTCGTCTGCGTCTCCAAAACAGAGAAGCCAACGCCAACGCCTCCATTTCGTCGAATCCCCTCCTCTCTT 
CTTCCATAATGTCAATTATCTC 
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GAGATAA.TTGnCATTATGGAAGAAGAGAGGAGGGGATTCGACGAAATGGAGGCGTTGGCGTTGGCTTCTCTGTTTT 

gagataattgHcattatggaagaagagaggaggggattcgacgaaatggaggcgttggcgttggcttctctgtttt 
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16 0 
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CGACAGACAAA 
CGACAGACAAA 
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CAGACGCCGTCTG 



CTTAGGGCTAC 
CTTAGGGCTAC 



D Open reading frame 



Mul3 
2mWD4 0 
SbWD4 0 
O3WD4 0 IJrC 



SDPNLQ 
SDPNLQ 
SDPNLQ 
SDPNLQ 



SSSATGGASP 
SSSATGGASP 
SSSATGGASP 
SSSATGGASP 



PSSFSGBCS 
PSS FSGECS 
PSS FSGBCS 
PSSFSGBCS 



Open reading frame 2 



Mul3 GPr 

ZmWD4 0 GpRt^ 

SbWD4 0 - -QtE 

OSWD40 GA 



EEGHIYSL 
EEGHIYSL 



Figure 1 Sequence characteristics of Mu13 element. (A) Sequences of Mu13. Terminal-inverted repeat (TIR) region is underlined. Bold sequences 
indicated two open reading frames and conceptual translation. (B) Structure and restriction map of Mu13. (C) Alignment of Mu13 left and right 
TIR. (D) Alignment of two ORFs of Mu13 element with three WD40 proteins, maize ZmWD40 (ACG25371 ), sorghum SbWD40 (Sb01 g008680), and 
rice OsWD40 (Os03g0738700). 



of Mul3y nine inbred lines of maize were analyzed by Southern blot 
analysis by using the Mul3 internal sequence as a probe. Six inbred 
lines (W22, M14, Q66, Q67, Q77, and Q79) that founded the Rob- 
ertson's Mutator population were included. To ensure a complete di- 
gestion, EcoRI was used, because it is methylation insensitive and does 
not cut Mul3 internally (Figure IB). As shown in Figure 3A, Mul3 
was detected in W22, Q67, Q77, and Q79, and was not detected in 
B73, Mol7, or A188. Mul3 was probably not present in M14 and Q66 
because the hybridized bands were substantially weak in comparison 
to other lines. PGR analysis by using Mul3 specific primers did not 
detect Mul3 in B73, Mol7, A188, M14, and Q66, suggesting that the 
weak signal may have arisen from non-specific hybridization with the 
probe. W22 appeared to contain two copies of Mul3, whereas other 



inbred lines contained one to two copies. In B73, a ~4.4kb fragment 
was weakly hybridized. This fragment is consistent with a WD40 gene 
(Accession no. ACG25371), which predicts a 4382 kb EcoRI fragment. 
The fragment contained a 258bp region that has 95% identity, and 
a 126bp region that has 98% identity to the Mul3 probe. It was 
expected to weakly hybridize with the Mul3 probe. In Mol7, the 
corresponding fragment is predicted to be 4302bp which was cross- 
hybridized as well. The Mul3 probe carried an 80bp sequence of the 
TIR (refer to File SI), which predictably would weakly hybridize with 
related Mu elements. This might explain the background and weak 
signals. 

To test whether Mul3 is present in the ancestor of maize, we 
analyzed Zea mays subsp. parviglumis (Accession no. PI 384061) 
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Figure 2 Detection of new transpositions of the Mu13 element in the UniformMu population. Genomic DMAs from eighteen randomly selected 
lines that segregate different seed mutant phenotypes were digested with EcoRI and hybridized with a Mu13 (A) and a Mu1 specific probe (B). A 
pooled WT (non-segregant, N) and a segregant (segregating each mutant phenotype, S) sample were used from each line (refer to Materials and 
Methods). Arrows indicate Mu13 insertions that were not found in the progenitors. 



and Zea mays subsp. mexicana (Accession no. PI 566684). The former 
is believed to be the ancestor of maize from a single domestication 
process (Matsuoka et al 2002). To reduce the chance that Mul3 may 
have resulted from a large fragment that escaped from Southern de- 
tection, five restriction enzymes {EcoRIy EcoRVy Hindllly Kpnl and 
Sail) that did not digest inside Mul3 were used to restrict the genomic 
DNA. Zea mays subsp. parviglumis did not contain any Mul3 ele- 
ment, as indicated by the absence of a hybridized signal (Figure 3B). A 
5kb Hindlll fragment was detected in Zea mays subsp. mexicanay but 
the signal intensity was much weaker than the Mul3 signals in W22 
or Q79. Because this hybridization was carried out under the same 



conditions at which the inbred DNAs were hybridized (and the load- 
ing was comparable with samples such as W22 or Q79), the signal 
produced in Zea mays subsp. mexicana was more likely from the 
WD40 fragment or an unknown homologous fragment than the real 
Mul3 element. Although the primers were proven robust, subsequent 
PGR detection by Mwi 3-specific primers failed to amplify the Mul3 
element from Zea mays subsp. mexicanay indicating that Mul3 was 
not present in Zea mays subsp. mexicana either. This result indicated 
that Mul3 is not present in the sample of two teosinte accessions 
tested, but as substantial genetic diversity exists among teosinte acces- 
sions, we cannot infer its absence among all teosintes. 
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Figure 3 Southern analysis of Mu13 
element presence in teosinte and 
maize lines. (A) Genomic DNAs of se- 
lected maize inbreds were digested 
with EcoRI and hybridized with a 
Mu13 probe. The arrow indicates the 
4382bp EcoRI fragment of the WD40 
gene (Accession no. ACG25371). (B) 
Genomic DNAs from Zea mays subsp. 
parviglumis (Accession no. PI 384061) 
and Zea mays subsp. mexicana (Acces- 
sion no. PI 566684) were digested with 
five different restriction enzymes (RV: 
EcoRV, H3: Hindlll, Kpnl, Sail, Rl: 
EcoRI) and hybridized with a Mu13 
probe. 
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Insertions of Mu13 in functional genes 

A Mul3 insertion was first identified in molecular characterization of 
the ZmSigS gene (Accession no. CG893004, GRMZM5G830932). The 
ZmSig3 gene consists of six exons, and the Mul3 was inserted in the 
third exon (Figure 4B). Analysis of the progenitor lines and a popula- 
tion segregating zmsigS mutants by PGR using Mu-TIR primer (TIR8) 
and the ZmSig3 specific primer (ZmSig3-R) indicated that this inser- 
tion was not present in WT and the progenitor lines, suggesting that it 
was a new transposition event. 

The insertion in ZmSigS suggested that the Mul3 element may 
contribute significantly to mutagenesis in the UniformMu population. 
In a previous study on seed mutants isolated from the UniformMu 
mutagenesis, Mu-flanking sequences were extracted by the Mu-TAIL 
method (deposited in GenBank, McCarty et al. 2005). The identity of 
the Mu element, however, is unknown. To search for insertions by 
Mul3, we chose two Mu insertions in known functional genes for 
analysis. Analysis of the Mu flanking sequences indicated that one 
Mu element was inserted in a paralog of Vpl4 gene on chromosome 
5S, named as Vpl4b (GRMZM5G838285), and that the other was 
inserted in a gene coding for an NAG (NAM, ATAF1,2, GUG2) 
transcription factor domain containing protein, named ZmNAGl 
(GRMZM2G3 12201). Both insertions were novel, as they were not 
present in the progenitors and were segregated specifically in individ- 
ual lines. We cloned and sequenced the inserted Mu elements. The 
insertion in Vpl4b was a Mu7 element (data not shown), and the one 
in ZmNAGl was a Mul3 element (Figure 4D). We analyzed a pop- 
ulation of 12 individual plants derived from a selfed heterozygote of 
the Mul3 insertion by using Mu TIR specific primer TIR8 and 
ZmNACl specific primers. The 12 individual plants were genotyped 
(Figure 4G). Plants homozygous for the Mul3 insertion showed 
a dwarf phenotype (#5, #8, and #10 in Figure 4E). Although the in- 
sertion was not confirmed as the cause of the dwarf phenotype, it 



indicated at least a linkage between this Mul3 insertion and the dwarf 
phenotype. The ZmNACl is likely a functional gene, as multiple ESTs 
were found in GenBank. 

Presence and transposition of six new Mu elements in 
UniformMu lines 

The sequencing of the maize genome revealed a surprising view of the 
Mutator family, which accounts for 1% of the B73 genome (Schnable 
et al 2009). We used the conserved 200bp Mu TIR sequences and 
performed a BLAST search of the maize genomic sequences in 
GenBank. A high stringency search (E value < e~^°) resulted in 
a low return of Mu elements. It appears that four types of Mu 
elements with distinct TIRs are present in the B73 genome. One 
class of Mu elements possesses TIRs with high similarity to known 
Mu elements in sequence and length (left and right TIR ~210 bp). 
A second class contains a left TIR of ~215bp and a short right TIR 
of ~90bp. A third class contains both short TIRs (~100bp), and 
a fourth class is called SOLOs, which contain only one TIR. The 
previously known Mu elements only account for a very small frac- 
tion of this family. Mwi, Mw2, Mw8, and Mul3 do not exist in the 
B73 genome. B73, however, does contain truncated and apparently 
non-functional derivatives of the autonomous MuDR, as well as 
one copy of Mu3 and Mw7, two copies of Mu4 and Mw5, and four 
copies of Mu7 derivatives that have insertions or deletions in their 
internal sequences. Because of the absence of MuDR, these ele- 
ments are dormant, and some may have lost their transposition 
activity due to accumulated mutations. To analyze the presence 
and possible activity of the unknown Mu elements in the 
UniformMu population, we identified a subset of MULEs from 
the B73 genome. The criteria are that the element 1) contains a 
highly conserved ~220bp TIR on both ends (>85% identity to 
consensus Mu TIR sequence); 2) contains perfect TIR ends 
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Figure 4 Insertion of Mu73 element 
into a gene coding for putative plastid 
sigma factor 3 [ZmSlgS) and an NAC 
domain-containing protein (ZmNACl). 
(A) PCR segregation analysis of 
] ZmSlgS mutant segregation popula- 
tion carrying a Mu73 element inser- 
tion. 1-17 were individual plants of 
an F2 family. The primers were TIR8- 
specific for Mu ends and ZmSig3-R- 
specific for ZmSlgS. (B) Structure of 
the ZmSlgS gene. Exons were desig- 
nated as boxes and introns as lines. 
The Mu73 insertion and primer anchor 
sites are indicated. The Mu^S element 
is not drawn to scale. (C) PCR genotyp- 
ing of a segregation population of the 
Mu1S insertion in ZmNACl. 1 -1 2 were 
individual plants of the population. 
ZmNACl -F and ZmNACl -R are pri- 
mers specific to ZmNACl. ++, +-, 
and — — designated WT, heterozy- 
gote and homozygote for the insertion 
in ZmNACl gene, respectively. The 
red arrow indicates a fragment con- 
taining a Mu73 element. (D) Structure 
of the ZmNACl gene and primer loca- 
tions. (E) Plants of #5, #8, and #10 (in 
C) exhibited a dwarf plant phenotype. 
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(GAGATA at the 5' and TATCTC at 3'); and 3) possesses perfect 
TSD in the insertion site, which is indicative for recent transposi- 
tion. The known active Mu elements all contain these features. We 
chose six MULEs that showed the highest similarity to Mu TIR 
consensus and with unrelated internal sequences. These elements 
were named as Mul4 to Mul9 (Accessions no. HQ698273- 
HQ698278, refer to File SI). 

A phylogenetic analysis performed by using the internal sequences 
indicated that these elements are not related, except for Mul and Mu2 
(also known as Mul.7), Mu5, and MuDR (Figure 5A). Sequence anal- 
ysis strongly supports the notion that Mul is a deletion derivative of 
Mw2, and that Mu5 a deletion derivative of MuDR. Mu2 contains 
a 140bp direct repeat in the internal sequence (Figure 5B, box arrows). 
Its 3' region (from 893 to 1330bp) showed a 91% identity to maize 
and rice genomic sequences, suggesting possible gene fragment cap- 
ture. The conceptual translation product of this region showed high 
similarity (80% identity) to OsOSgO 128200, which was annotated as 
zinc finger CCCH domain- containing protein 33 in rice. In sequence 
alignment with Mw2, Mul lacks most of this region, but still retains 



a residual 41 bp of the likely captured fragment. Similarly, Mu5 con- 
tains two segments of the mudrA gene that codes for MURA trans- 
posase (Figure 5B). Mul5 showed slight similarity to Mul9, in which 
three short segments of the internal sequences shared some similarity, 
suggesting that the two elements are likely of the same origin. The 
divergent sequences indicated that deletion and insertion also oc- 
curred fairly long ago. It has been known that the internal sequences 
of Mu elements are likely captured gene fragments. The captured gene 
fragments were analyzed in the known Mu elements (Lisch 2002). Our 
analysis indicated that Mu3 and Mu4 can be classified as Pack- 
MULEs. A fragment from maize chromosome 6 accounted for most 
of the internal sequence of Mu3. Additionally, two fragments fused 
from maize chromosome 1 and 3 accounted for the internal sequence 
of Mu4. 

We analyzed the internal sequences of new Mu elements identified 
in this work. As indicated in Figure 5B, the Mul 3 element contains 
two regions that showed high similarity to a WD40 protein. Mul4 
contains a fragment highly similar (89% identity) to a putative 
cucumisin-like serine protease on chromosome 1. Mul 5 contains 
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Figure 5 Sequence and structure of 
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a fragment that is similar to xylem serine proteinase 1 on chromo- 
some 1 (LOC100281759). Mul6 contains a fragment similar to a re- 
ceptor protein kinase TMK precursor (95% identity, Accession no. 
BT054484) on chromosome 3. Mul7 contains fragments from dif- 
ferent chromosomes. Mul8 contains a fragment of an auxin re- 
sponse factor 15 (ARF15) gene (Accession no. HM004530, 97% 
identity at nt level) and a calmodulin (LOC100286292, 98% identity 
at nt level). Mul9 contains a fragment of a putative xylem serine 
proteinase 1 (Accession no. NM_001154679). All of these elements 
except MuDR are between 1.4 and 2.5 kb in length. 

We cloned the internal sequences of these Mu elements and used 
them as probes to test their presence in inbred lines Mo 17, W22, and 
six randomly selected UniformMu mutant lines. As shown in a South- 
ern blot analysis (Figure 6), Mul4 to Mul8 elements were found in 
W22 and the UniformMu lines. The identical sizes of the fragments 
between W22 and the UniformMu lines strongly suggested that these 
elements were likely derived from W22. Mul9 was not found in 
Mol7, W22, or the UniformMu lines, but was found in B73. The 
analysis revealed that these elements represent part of the non-co- 
linear genome fraction of the three inbred lines. B73, Mo 17, and W22 
were all variable for these six elements in terms of copy numbers and 
RFLP size. Mul9 was not present in either Mo 17 or W22. It was also 
not detected in the six UniformMu lines. Some elements showed 
identical size among the three inbred lines, indicating likely early 
transposition events prior the separation of these inbreds. These Mu 
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elements in the UniformMu population were derived from W22. 
Within the limited number of the UniformMu samples, new trans- 
position events were not detected. 

DISCUSSION 

Mu13 is a new Mu transposable element 

Mutator elements share a highly conserved ~220 bp TIR sequence 
and create a 9bp TSD upon insertion (Walbot and Rudenko 2002). 
Different Mu elements are defined by the internal sequences between 
the TIRs. Mul3 has a TIR of 223bp that is highly similar to the 
conserved TIR sequences of known Mu elements (Figure 1), and yet 
the internal sequence is completely different from known Mu ele- 
ments. Mul3 element was not found in the sequenced B73 genome, 
nor was it detected by Southern hybridization analysis (Figure 3). Of 
the two Mul3 insertions identified in this study, each created a 9bp 
TSD. Hence, we concluded that Mul3 is a new Mu element. 

Mu13 contributes significantly to mutagenesis in the 
UniformMu population 

Among the previously identified Mu elements, not all are equally 
active in transposition. Mu element transposition was driven by the 
autonomous element MuDR (Hershberger et al 1991). Mu4 and Mu5 
were found inactive (Talbert et al. 1989), which may likely be due to 
the absence of the MuDR element. However, in a large scale tagging of 
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Figure 6 Presence of the new 
Mu elements in W22, Mo17, 
and the UniformMu line. Geno- 
mic DNAs were digested with 
EcoRl and probed with the 
internal sequence of each Mu 
element (as indicated under- 
neath). UniM-1 to 5 were ran- 
domly selected mutant lines 
from the UniformMu popula- 
tion. Fragment sizes are indi- 
cated by arrows. 
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gl8 locus, 58% insertions were caused by MullMu2, 25% by MuDR, 
7% by Mull, and the remaining elements {Mu8, Mul2, Mu4 and 
MulO) collectively merely contributed to 10% (Dietrich et al 2002). 
Our analysis of the known Mu element transposition events in the 
UniformMu population showed somewhat different presentation. 
Mul showed higher frequency of transposition, followed by Mu3, 
Mu8, and MuDR (B. C. Tan and D. R. McCarty, unpublished data). 

We have provided evidence that Mul 3 is active by analyzing new 
transposition events in a random selection of the UniformMu mutant 
lines, as well as by identification of Mul 3 insertions in two functional 
genes. After the initial identification of Mul 3, we recovered another 
Mul 3 insertion in one of two insertions analyzed. These results in- 
dicated that Mul 3 is highly active in transposition in the UniformMu 
population. Conceivably, if Mul 3 is present in other Mu active lines, it 
should contribute significantly to mutagenesis as well. All the Mu 
active lines were derived from a single line from which Mu transpos- 
able elements were discovered (Chandler and Hardeman 1992). The 
UniformMu population was derived by introgressing Robertson's Mu 
active line into inbred W22. Robertson maintained the Mu activity by 
out-crossing with W23; thereby Mul3 could have been derived from 
W23. However, our Southern blot analysis as well as PCR detection 
confirmed the presence of Mul 3 in W22. Comparison of the fragment 
sizes of Mul 3 inserted elements between W22 and the UniformMu 
lines indicated that they are identical. Unless W23 contained the 
similar Mul3 insertions, the Mul3 elements in the UniformMu pop- 
ulation appeared to be derived from W22. Hence, it is possible that the 
introgression of MuDR elements from the Mu active line activated the 
Mul 3 element in the W22. It would be interesting to test the presence 
of the Mul 3 element in other mutagenesis populations, such as the 
Maize-targeted mutagenesis population (May et al. 2003) and the 
Pioneer Hi- Bred International's Trait Utility System in Corn collec- 
tion (Bensen et al 1995). 

Unidentified Mu elements in the maize genome and 
their activity 

The identification of active Mul 3 in the maize genome suggested that 
there are many more unknown Mu elements in the genome. The 
sequencing of the B73 genome recovered many of these elements 
(Schnable et al 2009), but evidence suggested that there are more. 
We have analyzed six Mu elements identified in B73, and Mul 9 was 
not detected in the UniformMu population, which is largely W22 
introgressed with the Robertson's active Mu line. Mul 3 was not found 
in the B73 genome, and in the tagging of 80 alleles of glS, Mul 3 was 
not detected in the population that derived from the Robertson's Mu- 
active line. Mul3 apparently was present in W22 and was activated 
during introgressing with Mu-active lines. The W22 line did not 
contain any active MuDR element (B. C. Tan and D. R. McCarty, 
unpublished data), hence Mul 3 was inactive. Because different maize 
inbred lines harbor different spectrums of Mu elements, more un- 
known Mu elements are expected. It is highly likely that most of the 
maize inbred lines did not contain any active MuDR elements, hence 
all the MuDR driven Mu elements are dormant. Upon introducing the 
MuDR element, Mu element activity may be restored. If this is the 
case. Southern blot based cosegregation analysis using known Mu 
internal sequences as probes may encounter some problems. But this 
will not affect analysis based on the TIR sequences such as AIMS 
(Frey et al 1998), Mu-TAIL PCR (Settles et al 2004), AIMS and 
Mu-TAIL-PCR combined (Yi et al 2009), and the use of PCR-coupled 
with pyrosequencing (Williams- Carrier et al 2010). In addition, if 
the creation of new Mu elements is associated with the MuDR activity, 
it will be expected that there will be many new Mu elements in 



Robertson's Mu-active line. It will be interesting to know the Mu 
landscape in the Robertson Mu active line. 
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